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ABSTRACT:  Starting  from  previous  benchmark  CBS-QB3  electronic  structure 
calculations  (Conforti,  P.  F.;  Braunstein,  M.;  Dodd,  J.  A.  /.  Phys.  Chem.  A  2009,  113, 
13752),  we  develop  two  global  potential  energy  surfaces  for  0(3P)  +  DMMP  collisions, 
using  the  specific  reaction  parameters  approach.  Each  surface  is  simultaneously  fit  along 
the  three  major  reaction  pathways:  hydrogen  abstraction,  hydrogen  elimination,  and 
methyl  elimination.  We  then  use  these  surfaces  in  classical  dynamics  simulations  and 
compute  reactive  cross  sections  from  4  to  10  km  s-1  collision  velocity.  We  examine  the 
energy  disposal  and  angular  distributions  of  the  reactive  and  nonreactive  products.  We 
find  that  for  reactive  collisions,  an  unusually  large  amount  of  the  initial  collision  energy 
is  transformed  into  internal  energy.  We  analyze  the  nonreactive  and  reactive  product 
internal  energy  distributions,  many  of  which  fit  Boltzmann  temperatures  up  to  ~2000 
K. 
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I.  INTRODUCTION 

The  proliferation  of  weapons  of  mass  destruction,  including 
chemical  warfare  agents,  has  renewed  interest  in  the 
fundamental  chemistry  of  nerve  agents  such  as  sarin  and  VX 
and  their  main  simulant,  dimethyl  methylphosphonate 
(DMMP).1-8  For  example,  recent  studies  have  examined  the 
interaction  of  DMMP  with  surfaces,  including  reactivity  and 
decomposition  on  Ti021  and  silica.2  There  have  also  been 
recent  investigations  on  gas  phase  decomposition  mechanisms 
and  reaction  rate  kinetics  of  DMMP.3,8'9  In  the  present  work, 
we  focus  on  the  interaction  of  DMMP  with  atomic  oxygen, 
0(3P),  which  is  important  in  assessing  chemical  interactions  of 
nerve  agents  released  in  the  upper  atmosphere,  where  0(3P)  is 
a  major  species. 

In  our  initial  study  of  0(3P)  collisions  with  sarin  and  the 
simulant  DMMP,5  we  characterized  the  main  stationary  points 
and  reaction  pathways  with  several  levels  of  theory,  including 
semiempirical,  density  functional,  and  benchmark  CBS-QB3. 
For  these  systems,  we  found  three  major  pathways.  For  0(3P) 
+  DMMP,  the  focus  of  the  present  work,  the  three  pathways  are 
shown  in  Figure  1:  hydrogen  abstraction  to  create  OH 
(reactions  1  and  2),  the  elimination  of  hydrogen  following 
0(3P)  addition  (reactions  3  and  4),  and  the  elimination  of  a 
methyl  radical  following  0(3P)  addition  (reactions  5  and  6). 
The  hydrogen  abstraction  reactions  have  relatively  low  reaction 
barriers  less  than  10  kcal  mol-1  and  reaction  enthalpies  that  are 
within  5  kcal  mol-1  of  thermoneutral.  In  contrast,  the 
elimination  reactions  have  barriers  greater  than  40  kcal  mol-1 
and  reaction  enthalpies  that  range  from  approximately  —15  to 
+35  kcal  mol-1.  These  characteristics  have  a  great  deal  in 
common  with  0(3P)  plus  small-chain  hydrocarbon  systems, 
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which  also  have  comparable  reactive  channels,  transition  states, 
and  reaction  energies.10-12 

Starting  from  our  previous  work  on  the  0(3P)  +  DMMP  and 
0(3P)  +  sarin  systems,  in  the  present  paper  we  develop  two 
global  potential  energy  surfaces  for  0(3P)  +  DMMP  collisions, 
using  the  specific  reaction  parameters  (SRP)  approach.13 
Truhlar  and  co-workers  initially  developed  the  SRP  method, 
where  semiempirical  potential  parameters  were  tuned  in  order 
to  examine  a  reaction  of  interest,  and  then  direct  dynamics  were 
performed  using  the  modified  semiempirical  parameters.  As 
many  groups  have  done  since  the  seminal  work  of  Truhlar,  we 
adopt  the  SRP  approach  and  use  it  to  examine  the  0(3P)  + 
DMMP  potential  surface.  Each  of  the  global  surfaces  is 
simultaneously  fit  using  benchmark  CBS-Q_ electronic  structure 
calculations  along  the  three  major  reaction  pathways.  We  then 
use  these  surfaces  in  classical  dynamics  simulations.  We 
compute  the  cross  sections  for  the  major  0(3P)  +  DMMP 
reaction  pathways,  from  4  to  10  km  s-1  collision  velocity.  We 
then  examine  the  energy  disposal  and  angular  distributions  of 
the  reactive  and  nonreactive  products.  We  find  that  for  reactive 
collisions,  an  unusually  large  amount  of  the  initial  collision 
energy  is  transformed  into  internal  energy.  For  example,  for  H- 
atom  elimination  at  8  km  s-1  collision  velocity,  on  average 
about  80%  of  the  available  energy  goes  into  internal  energy  of 
the  product  DMMP— O  radical.  This  ^percentage  is  larger  than 
comparable  hyperthermal  reactions14  and  many  chemical 
interactions  in  general.  For  CH3  elimination,  about  60—70% 
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Figure  1.  Major  collisional  reactions  for  DMMP  +  0(3P):  (l)  and  (2) 
hydrogen  abstraction;  (3)  and  (4)  hydrogen  elimination;  (5)  and  (6) 
methyl  elimination. 


goes  into  internal  energy  of  products,  and  for  H-atom 
abstraction  between  50  and  60%  goes  into  product  internal 
energy.  The  internal  energy  deposition  is  correlated  with  the 
distribution  of  impact  parameters:  nearly  all  H-elimination 
reactions  are  “hard  collisions”  with  impact  parameters  less  than 
1.0  A.  Nearly  all  CH3  elimination  reactions  occur  for  impact 
parameters  less  than  2  A,  and  the  H-atom  abstraction  pathway 
has  a  broad  range  of  impact  parameters  up  to  about  3  A. 
Nonreactive  collisions  result  in  a  wide  range  of  product  internal 
energies,  which  can  be  separated  into  forward  and  backward 
scattering  mechanisms.  We  further  analyze  the  nonreactive  and 
reactive  product  internal  energy  distributions,  many  of  which  fit 
Boltzmann  temperatures  up  to  ~2000  K,  and  predict  excited 
product  spectra. 

Hyperthermal  O-atom  beams  have  recently  been  used  to 
investigate  0(3P)  +  small-chain  hydrocarbon  reactions,10,11,15,16 
and  there  have  been  recent  studies  which  use  SRP  potentials  to 
examine  these  same  systems.12,17  As  we  are  not  aware  of  any 
other  efforts  examining  hyperthermal  collisions  of  0(3P)  with 
DMMP  or  other  nerve  agents  or  simulants,  the  present  results 
provide  guidance  for  possible  future  measurements  in  this  area. 
Furthermore,  the  present  results  extend  the  use  of  SRP 
potentials  to  larger  hydrocarbon  systems  and  give  insight  into 
the  unusual  hyperthermal  collision  dynamics  of  medium-sized 
hydrocarbon  molecules. 

The  paper  proceeds  as  follows.  Section  II  describes 
benchmark  electronic  structure  calculations  performed  on  the 
0(3P)  +  DMMP  system  and  the  construction  of  two  separate 
global  potential  energy  surfaces  using  the  SRP  potential  form. 
The  construction  of  two  potential  surfaces,  which  use  slightly 
different  cost  functions,  is  useful  as  it  gives  some  indication  of 
the  sensitivity  of  the  cross  sections  and  other  dynamical 
quantities  to  errors  in  the  surface.  The  surfaces  are  compared  to 
semiempirical  potentials  and  to  the  benchmark  calculations 


from  which  they  were  derived.  Section  III  describes  0(3P)  + 
DMMP  classical  dynamics  calculations  using  the  two  SRP 
surfaces.  Reaction  cross  sections  are  given,  and  product  internal 
energies  and  angular  distributions  are  analyzed.  Section  IV  gives 
conclusions. 

II.  AMI  SRP  POTENTIAL 

Although  there  have  been  great  strides  in  developing  accurate, 
multidimensional,  global  potential  energy  surfaces,1  -24  mole¬ 
cules  with  more  than  approximately  six  atoms  remain  a 
challenge,  principally  due  to  the  large  number  of  degrees  of 
freedom  involved  and  the  possibility  of  multiple  reaction 
pathways.  Recently  the  SRP  method  has  been  developed  and 
successfully  applied  to  hyperthermal  reaction  systems  contain¬ 
ing  five  to  nine  atoms.12,17  The  idea  is  to  start  with  a 
semiempirical  potential  and  adjust  the  parameters  to  match 
properties  of  known  experimental  values  and/or  high-order 
electronic  structure  calculations  for  a  particular  process  or 
reaction.  The  fast-running,  full  degree  of  freedom  properties  of 
the  semiempirical  method  are  retained,  but  the  SRP  surface 
provides  accuracy  comparable  to  the  potential  surface  data  used 
in  the  reparameterization.  The  challenge  of  this  method  is  to 
optimize  the  semiempirical  parameters  in  such  a  way  that  key 
features  of  the  surface  are  improved  to  a  satisfactory  degree, 
without  introducing  deficiencies  and  artifacts. 

We  seek  to  optimize  a  single  semiempirical  potential  which 
simultaneously  describes  all  three  pathways  (six  reactions  of 
Figure  l)  for  0(3P)  +  DMMP  collisions.  We  chose  a  procedure 
that  has  been  used  for  other  smaller  systems,  including  alkane 
rearrangement  reactions,25,26  hyperthermal  0(3P)  +  alkane 
reactions,12,17  and  hyperthermal  F  +  alkane  reactions.27  We 
start  with  an  AMI  potential.  We  interfaced  the  nonlinear  least- 
squares  solver  MINPACK28  with  the  GAMESS  computational 
chemistry  package.29  While  other  algorithms  for  nonlinear 
optimization  of  semiempirical  parameters  have  been  utilized, 
such  as  genetic  algorithms,1  '  the  Levenberg— Marquardt 
algorithm  as  implemented  in  MINPACK  also  performs  well  in 
optimizing  parameters.17,31  For  optimization  of  the  AMI 
parameters,  we  chose  a  set  of  relevant  structures.  In  particular, 
the  reactants,  products,  and  transition  states  as  shown  in  Figure 
1  of  ref  5,  and  their  energies  and  gradients  at  the  CBS-QB3 
level  are  used.  The  combustion  study  of  Tester  and  co-workers 
showed  that  the  differences  between  the  ab  initio  calculations 
using  CBS-Q^of  the  heat  of  formation  for  organophosphorus 
compounds  were  within  2.5  kcal  mol-1  of  the  literature  values 
where  available.32  Therefore,  we  use  this  level  of  theory  as  the 
benchmark  for  the  optimization  calculations.  Additionally, 
starting  from  the  transition  structures,  we  used  10  points 
(energies  and  gradients)  in  both  the  forward  and  reverse 
directions  along  the  intrinsic  reaction  coordinate  using  B3LYP/ 
CBSB7,  which  is  the  initial  calculation  for  the  compound  CBS- 
0133  method  used  in  ref  5.  GAMESS  performs  the  semi¬ 
empirical  calculations  and  MINPACK  adjusts  AMI  parameters 
on  the  relevant  atoms  (P,  C,  H,  and  O)  using  the  least-squares 
cost  function 

V  / 17 AMI  cCBS\2  .  V  (t  AMI  t  CBS\2 

WE  2j  ( Ei  ~  Ei  )  +wLZj  ( Li  ~  Li  ) 

i  i 

where  E{  and  L,  are  the  energy  and  gradient  magnitudes, 
respectively,  of  geometry  i  with  weighting  wE  and  wL.  For  these 
calculations,  wE  is  set  to  1.0  (mol  kcal-1)2  and  wL  is  set  at  50.0 
(mol  A  kcal-1)2.  For  phosphorus,  oxygen,  and  carbon,  12  AMI 
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Table  1.  AMI  and  AM1-SRP  Parameters  for  H,  C,  O,  and  P  Atoms 


H 

C 

parameter 

AMI 

AMI -SRP  1 

AM1-SRP2 

AMI 

AM1-SRP1 

AM1-SRP2 

Uss 

-11.3964 

-11.7606 

-11.7423 

-52.0287 

-51.5477 

-51.5311 

upp 

-39.6142 

-38.8044 

-38.7211 

A 

-6.1738 

-6.2808 

-6.2837 

-15.7158 

-16.0314 

-16.0471 

-7.7193 

-7.8639 

-7.8835 

Cs 

1.1881 

1.1725 

1.1704 

1.8087 

1.7612 

1.7596 

4 

1.6512 

1.6749 

1.6732 

a 

2.8823 

2.9116 

2.9129 

2.6483 

2.6402 

2.6402 

Gss 

12.8480 

12.8918 

12.9121 

12.2300 

12.8842 

12.8535 

GSP 

11.1400 

11.8874 

11.8968 

Gpp 

11.0800 

10.6964 

10.7772 

GP2 

9.8400 

10.1710 

10.1809 

Hs  p 

2.4300 

2.4988 

2.4854 

O 

P 

parameter 

AMI 

AM1-SRP1 

AM1-SRP2 

AMI 

AM1-SRP1 

AM1-SRP2 

17SS 

-97.8300 

-98.7177 

-98.5818 

-42.0299 

-42.4211 

-42.4248 

Upp 

-78.2624 

-78.2643 

-78.2622 

-34.0307 

-34.5494 

-34.5681 

Ps 

-29.2728 

-28.9311 

-28.8875 

-6.3538 

-6.3010 

-6.2943 

P, 

-29.2728 

-29.1341 

-29.1139 

-6.5907 

-6.2516 

-6.2376 

fs 

3.1081 

3.1540 

3.1592 

1.9813 

2.0817 

2.0883 

4 

2.5240 

2.5107 

2.5100 

1.8752 

1.9499 

1.9548 

a 

4.4537 

4.4333 

4.4450 

2.4553 

2.4650 

2.4656 

Gss 

15.4200 

14.9095 

14.8658 

11.5600 

11.4228 

11.4388 

Gsp 

14.4800 

14.3813 

14.3836 

5.2374 

5.0413 

5.0344 

Gpp 

14.5200 

14.5104 

14.5002 

7.8776 

7.1561 

7.1431 

GP2 

12.9800 

12.9777 

12.9780 

7.3076 

6.9030 

6.8879 

Hs  P 

3.9400 

3.9528 

3.9447 

0.7792 

0.7515 

0.7508 

parameters  are  optimized,  and  for  hydrogen  5  AMI  parameters 
are  optimized.  Therefore  there  are  41  ( 12* (phosphorus, 
oxygen,  and  carbon)  +  5  hydrogen)  AMI  parameters  fit  with 
a  total  of  320  data  points:  23  geometries*2(energy  +  gradient 
magnitude)  for  reactions  1—4  and  6,  plus  45  geometries*2- 
(energy  +  gradient  magnitude)  for  reaction  5.  To  understand 
the  sensitivity  of  the  resulting  SRP  potential  to  the  fitting 
method,  two  AM1-SRP  potentials  were  determined  using  the 
above  procedure  (labeled  AM1-SRP1  and  AM1-SRP2).  For  the 
AM1-SRP1  potential,  all  six  reactions  were  used  at  one  time  to 
reparameterize  the  AMI  parameters.  For  the  AM1-SRP2 
potential,  AMI  parameters  were  optimized  sequentially.  First 
the  abstraction  reactions  were  optimized,  followed  by  the 
hydrogen  elimination  reactions,  and  finally  the  methyl 
elimination  reactions.  During  the  optimization,  the  parameters 
were  allowed  to  vary  within  10%  of  their  original  values.  A 
comparison  of  the  AMI  parameters  is  given  in  Table  1  for  the 
H,  C,  O,  and  P  atoms.  The  average  percent  change  between  the 
SRP  and  original  AMI  parameters  is  approximately  2.2%,  with  a 
maximum  change  of  9.3%. 

Table  2  compares  the  energies  of  the  stationary  points  and 
transition  structures  for  the  0(3P)  +  DMMP  reactions  using 
CBS-QB3,  which  is  considered  benchmark,  and  the  unmodified 
AMI,  AMI -SRP  1,  and  AM1-SRP2  potentials.  All  energies  are 
referenced  to  the  0(3P)  +  DMMP  reactant  asymptote,  and  the 
units  for  all  energies  are  kcal  mol-1.  The  CBS-QB3  calculations 
are  performed  using  the  Gaussian03  computational  chemistry 
package,33  and  the  AMI  semiempirical  calculations  are 
performed  using  the  GAMESS  package.29  One  major  deficiency 
of  the  AMI  potential  is  that  the  reaction  energies  are  more 
exothermic  than  the  CBS-QB3  results,  particularly  for  the 


hydrogen  abstraction  reaction.  This  problem  is  largely 
corrected  with  the  SRP  potentials.  In  fact,  the  SRP2  potential 
for  reaction  2  reproduces  the  CBS-QB3  results  well  for  the 
barrier  and  reaction  energies.  However,  the  SRP2  improvement 
is  modest  at  best  for  hydrogen  elimination,  while  SRP1  does 
well  for  hydrogen  elimination  reaction  4.  For  methyl 
elimination,  the  SRP  potentials  improve  the  initial  barrier  for 
reaction  5  and  the  barrier  for  reaction  6,  but  for  some  of  the 
reaction  energies,  the  SRP  potentials  are  slightly  more 
exothermic  than  AMI.  This  degree  of  optimization  with  the 
SRP  potentials  is  similar  to  recent  studies  for  0(3P)  + 
hydrocarbons.12'17  For  instance,  in  Yan  et  al.,  two  semiempirical 
SRP  potential  sets  were  examined  for  the  0(3P)  +  ethane 
system,  where  either  6  or  12  reactions  on  the  energy  surface  are 
used  for  the  reparameterization.12  The  agreement  between  the 
SRP  results  and  the  high-level  MRCI  calculations  for  the 
activation  barriers  and  reaction  energies  span  from  approx¬ 
imately  0  to  16  kcal  mol-1.  This  shows  the  level  of  agreement 
for  a  diverse  reaction  set  and  is  similar  to  what  is  achieved  with 
the  present  study.  As  will  be  shown,  the  dynamical  results  for 
SRP1  and  SRP2  generally  agree  well  with  each  other  and  are 
much  less  reactive  than  the  AMI  results,  especially  near 
thresholds.  This  indicates  that  the  optimization  procedures 
have  qualitatively  changed  the  AMI  potential  and  that 
dynamically  the  SRP1  and  SRP2  potentials  are  behaving  in  a 
similar  way. 

To  compare  the  potential  surfaces,  we  calculate  the  absolute 
error  for  the  semiempirical  surfaces  using  discrete  points  along 
the  reaction  pathway  for  each  reaction  compared  to  the 
B3LYP/CBSB7  reaction  surfaces.  These  points  are  computed 
using  the  transition  state  geometry  and  scanning  those  bonds 
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that  are  breaking/forming  for  a  particular  reaction  while  leaving 
the  remaining  degrees  of  freedom  at  their  initial  geometry. 
Histograms  of  the  absolute  errors  for  points  along  the 
minimum  energy  for  all  reactions  are  given  in  Figure  2a  for 
the  nonoptimized  AMI  (black),  the  AM1-SRP1  (red),  and  the 
AM1-SRP2  (blue).  From  this  plot,  the  original  AMI  potential 
leads  to  frequent  errors  of  approximately  12.5—17.5  kcal  mol-1. 
Many  of  these  errors  are  in  the  product  channels  where  the 
AMI  is  overly  exothermic.  Both  the  AM1-SRP1  and  AM1- 
SRP2  decrease  the  large  errors  associated  with  the  AMI; 
however,  the  AM  1-SRP2  does  have  several  large  errors  (greater 
than  20  kcal  mol-1).  The  overall  mean  absolute  error  (MAE) 
for  the  nonoptimized  AMI  is  8.6  kcal  mol-1,  while  the  MAEs 
for  SRP1  and  SRP2  are  8.0  and  8.2  kcal  mol-1,  respectively. 
The  best  improvement  occurs  for  the  hydrogen  abstraction 
reactions  which  dominate  the  dynamics  at  lower  velocities.  In 
Figure  2b,  histograms  are  given  for  the  hydrogen  abstraction 
reaction  for  the  nonoptimized  AMI  (black),  the  AM  1-SRP1 
(red),  and  the  AM1-SRP2  (blue).  Similarly,  the  larg  e  error 
between  12.5  and  17.5  kcal  mol-1  for  the  original  AMI 
potential  is  due  to  the  large  exothermic  inaccuracy  for  these 
reactions  (as  shown  in  Table  2).  Additionally,  the  AM  1-SRP2 
shows  the  best  improvement  along  the  hydrogen  abstraction 
reaction  pathways.  For  this  reaction,  the  MAE  decreases  from 
8.7  kcal  mol-1  for  the  original,  unmodified  AMI  to  5.8  and  3.6 
kcal  mol-1  for  the  SRP1  and  SRP2,  respectively.  To  get  more  of 
a  sense  of  the  errors  in  the  surfaces,  contour  plots  of  the 
hydrogen  abstraction  reaction  1  are  given  in  Figure  3  and  plots 
of  the  hydrogen  abstraction  reaction  2  are  given  in  Figure  4  for 
each  level  of  theory.  The  overall  correction  of  the  overly  large 
AMI  exothermicities  is  clearly  seen  in  the  reparameterized 
surfaces. 

III.  MOLECULAR  DYNAMICS  SIMULATIONS 

Classical  molecular  dynamics  simulations  of  0(3P)  +  DMMP 
collisions  are  performed  using  the  dynamic  reaction  path 
(DRC)  routine  within  the  GAMES S  computational  chemistry 
package.29  The  potential  is  calculated  as  the  simulation 
progresses  using  both  AM1-SRP1  and  AM1-SRP2  potentials. 
For  reference,  we  also  include  results  using  the  standard, 
unmodified  (non-SRP)  AMI  potential.  Procedures  follow  those 
of  ref  5.  The  integration  time  step  was  0.1  fs  using  the  velocity- 
verlet  integrator,  the  initial  separation  between  fragments  was 
15  A,  final  separation  to  products  was  18  A,  the  SCF 
convergence  criterion  was  1.0  X  10-5  hartree,  and  the 
maximum  tolerance  for  the  deviation  for  energy  conservation 
(per  step)  is  set  to  5.0  X  10-5  hartree.  Collision  velocities  of  4, 
6,  8,  and  10  km  s-1  are  examined.  This  velocity  range  overlaps 
with  recently  reported  hyperthermal  measurement  capabil¬ 
ities34  and  also  overlaps  the  relative  collision  velocity  of 
outgassed  molecules  from  spacecraft  and  the  atmospheric  wind 
at  altitudes  of  150—700  km.35  At  each  collision  velocity  and  for 
each  surface  (the  standard  AMI,  AM1-SRP1,  AM1-SRP2), 
10  000  trajectories  are  calculated.  For  a  given  collision  velocity, 
at  most  40  trajectories  per  10  000  do  not  meet  the  maximum 
energy  tolerance  and  are  discarded.  To  start  a  trajectory,  the 
minimum  energy  structures  of  DMMP  are  randomly  oriented 
using  Euler  angles  and  a  random  impact  parameter  is  chosen. 
Initial  internal  DMMP  velocities  are  determined  from  the  zero- 
point  energy  through  normal-mode  sampling  with  zero  initial 
rotational  energy.  The  maximum  impact  parameter  is  set  at  4.1 
A.  Examination  of  the  opacity  function  has  shown  this  to  be  an 
acceptable  cutoff  value.  Most  of  the  trajectories  are  run  for 
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Figure  2.  Histograms  of  the  absolute  errors  for  points  along  the  energy  reaction  pathways  for  (a)  all  reactions  and  (b)  hydrogen  abstraction 
reactions.  The  AMI  results  are  shown  in  black,  the  AM1-SRP1  results  are  shown  in  red,  and  the  AM1-SRP2  results  are  shown  in  blue.  Error  bars 
indicate  standard  deviation,  and  the  energy  bin  size  is  2.5  kcal  mol-1. 


roh/A  roh/A 


Figure  3.  Surface  contours  for  the  hydrogen  abstraction  reaction  1  using  (a)  AMI,  (b)  AM1-SRP1,  (c)  AM1-SRP2,  and  (d)  B3LYP/CBSB7.  The 
energies  are  given  in  relation  to  the  DMMP  +  0(3P)  reactant  asymptote. 


several  hundred  femtoseconds  to  one  picosecond  to  assess  the 
primary  collision  interaction  process  and  resulting  reaction 
products. 

The  total  cross  section  as  a  function  of  collision  velocity  is 
shown  in  Figure  5  for  each  pathway.  In  the  plot  and  in  the 
following  discussion,  the  two  reactions  for  each  pathway  are 
summed,  so  that  reactions  1  and  2  of  Figure  1  are  summed  to 
give  the  H-abstraction  cross  section,  reactions  3  and  4  give  the 
H-elimination  cross  section,  and  reactions  5  and  6  give  the  CH3 
elimination  cross  section.  Three  results  are  shown  using  the 
three  different  potentials,  AMI  which  we  reported  on 
previously,5  SRP1,  and  SRP2.  In  general,  both  AM1-SRP 
cross  sections  are  similar  to  each  other  and  are  lower  than  the 


original  AMI  cross  sections  previously  calculated,  particularly 
the  H-abstraction  and  CH3  elimination  pathways.  The  major 
correction  of  the  SRP  potentials  is  to  decrease  the  reaction 
exothermicity  and  to  generally  increase  the  potential  energy 
after  the  barrier  region  in  the  product  channels.  The  increase  in 
the  potential  in  this  way  makes  it  less  likely  that  a  given 
trajectory  will  proceed  to  separated  products.  Focusing  on  the 
individual  reactive  pathway  SRP  cross  sections,  the  H-atom 
abstraction  reaction  has  a  significant  cross  section  at  low 
collision  velocities,  as  might  be  expected  by  the  relatively  low 
reaction  barriers.  The  hydrogen  elimination  reaction  becomes 
competitive  with  the  abstraction  reaction  at  8  km  s-1  when  the 
available  collision  energy  (~108  kcal  mol-1)  is  well  above  the 
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energies  are  given  in  relation  to  the  DMMP  +  0(3P)  reactant  asymptote. 
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Figure  5.  Reaction  cross  section  as  a  function  of  collision  velocity  (km 
s-1)  and  energy  (kcal  mol-1)  for  the  hydrogen  abstraction  (black), 
hydrogen  elimination  (red),  and  methyl  elimination  (blue)  reactions. 
The  AMI  results  are  shown  with  the  solid  lines  and  square  symbols; 
the  AM  1-SRP1  results  are  shown  with  dashed  lines  and  circular 
symbols;  the  AM  1-SRP2  results  are  shown  with  dotted  lines  and 
triangular  symbols.  Error  bars  indicate  standard  deviation. 


reaction  barrier.  At  lower  collision  velocities,  the  hydrogen  and 
methyl  elimination  pathways  are  comparable.  The  H- 
elimination  cross  sections  grow  much  more  rapidly  with 
collision  velocity  than  CH3  elimination,  so  that  at  10  km  s_1,  H- 
elimination  is  about  10  times  more  likely  than  CH3  elimination. 


Some  of  this  difference  is  due  to  the  fact  that  there  are  three 
times  as  many  available  hydrogen  atoms  as  methyl  groups. 

The  distribution  of  energy  in  the  collided  0(3P)  +  DMMP 
products  has  not  been  previously  explored.  The  product 
energies  will  be  important  in  guiding  product  detection 
schemes  for  possible  hyperthermal  beam  measurements. 
These  distributions  are  also  of  fundamental  importance  in 
understanding  the  collision  dynamics  of  O-atom  with  medium¬ 
sized  hydrocarbon  molecules.  Tables  3—5  present  an  overview 


Table  3.  Average  Total  Vibrational  Energy  (£vib)  Excluding 
Zero-Point,  Rotational  ( Erot ),  and  Relative  Translational 
Energy  (Etrans),  in  kcal  mol-1  for  the  Products  of  the 
Hydrogen  Abstraction  Pathwaya 

H  abstraction 


(km  s  ) 

potential 

(Fvib) 

(Trot) 

(Ttrans) 

8 

AMI 

52.1  (0.35) 

40.8  (0.27) 

56.7  (0.38) 

AM1-SRP1 

47.5  (0.36) 

38.7  (0.29) 

48.2  (0.36) 

AM1-SRP2 

34.7  (0.27) 

32.7  (0.25) 

63.0  (0.48) 

10 

AMI 

59.8  (0.28) 

48.9  (0.23) 

104.7  (0.49) 

AM1-SRP1 

58.0  (0.28) 

51.0  (0.25) 

94.9  (0.46) 

AM1-SRP2 

45.8  (0.24) 

40.3  (0.21) 

107.8  (0.55) 

aThe  fraction  of  the  total  energy  (initial  translational  energy  +  the 
average  of  the  reaction  enthalpies)  is  given  in  parentheses. 


of  the  product  energy  deposition  for  the  hydrogen  abstraction, 
hydrogen  elimination,  and  methyl  elimination  pathways, 
respectively.  The  average  vibrational  (excluding  zero-point), 
rotational,  and  relative  translational  energies  for  each  pathway 
are  given  at  collision  velocities  of  8  and  10  km  s_1.  As  in  the 
discussion  of  the  cross  sections,  results  for  the  two  channels  of 
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Table  4.  Average  Total  Vibrational  Energy  (E^h)  Excluding 
Zero-Point,  Rotational  ( Erot ),  and  Relative  Translational 
Energy  ( Etrans ),  in  kcal  mol-1  for  the  Products  of  the 
Hydrogen  Elimination  Pathway*1 

H  elimination 


velocity 


(km  s  x) 

potential 

(Fvib) 

<Trot> 

(-^trails) 

8 

AMI 

78.1  (0.62) 

22.9  (0.18) 

23.9  (0.19) 

AMl-SRPl 

73.9  (0.67) 

17.1  (0.16) 

18.8  (0.17) 

AM1-SRP2 

84.0  (0.63) 

29.4  (0.22) 

20.7  (0.15) 

10 

AMI 

115.9  (0.59) 

47.6  (0.24) 

32.3  (0.17) 

AMl-SRPl 

115.4  (0.63) 

43.4  (0.24) 

23.3  (0.13) 

AM1-SRP2 

129.2  (0.65) 

44.6  (0.22) 

25.9  (0.13) 

aThe  fraction  of  the  total  energy  (initial  translational  energy  +  the 
average  of  the  reaction  enthalpies)  is  given  in  parentheses. 


Table  5.  Average  Total  Vibrational  Energy  (E^h)  Excluding 
Zero  Point,  Rotational  (Erot),  and  Relative  Translational 
Energy  (Etrans)>  *n  kcal  mol-1  for  the  Products  of  the  Methyl 
Elimination  Pathway** 

CH3  elimination 


(km  s~{) 

potential 

<ivb> 

(Trot) 

(Ttrans) 

8 

AMI 

66.3  (0.52) 

14.0  (0.11) 

48.1  (0.38) 

AMl-SRPl 

61.9  (0.49) 

14.4  (0.11) 

49.1  (0.39) 

AM1-SRP2 

60.4  (0.48) 

17.4  (0.14) 

47.5  (0.38) 

10 

AMI 

97.6  (0.55) 

17.9  (0.10) 

62.0  (0.35) 

AMl-SRPl 

94.1  (0.52) 

22.5  (0.12) 

64.1  (0.35) 

AM1-SRP2 

109.6  (0.58) 

26.7  (0.14) 

50.1  (0.27) 

aThe  fraction  of  the  total  energy  (initial  translational  energy  +  the 
average  of  the  reaction  enthalpies)  is  given  in  parentheses. 


each  reaction  pathway  are  averaged  together.  In  Tables  3—5,  we 
also  give  the  fraction  of  total  energy  deposited  in  vibration, 
rotation,  and  translation  in  parentheses.  The  total  energy  is 
defined  by  the  sum  of  the  collision  energy  and  the  average 
enthalpy  of  the  reaction.  For  H-atom  elimination,  an  unusually 
large  amount  of  the  total  energy  (about  80%)  goes  into  internal 
energy  of  the  product  DMMP— O  radical.  For  CH3  elimination 
60—70%  goes  into  product  internal  energy,  and  for  H-atom 
abstraction  between  50  and  60%  goes  into  product  internal 
energy.  For  the  elimination  reactions,  most  of  this  internal 
excitation  goes  into  vibration.  As  previously  observed  in  the 
literature,12'17  H-abstraction  at  these  velocities  occurs  by  a 
“stripping”  mechanism  involving  glancing  collisions.  H- 
abstraction  therefore  does  not  redistribute  much  of  the 
available  energy  into  internal  excitation,  whereas  H-elimination 
and  to  a  lesser  extent  CH3  elimination  involve  low-impact 
parameter,  or  “hard”,  collisions.  Through  an  examination  of  the 
opacity  function  P(b),  these  mechanisms  are  readily  apparent. 
The  opacity  function  has  the  following  relation  to  the  total 
reactive  cross  section,  a  -  2nJ^m3X  P(b)b  d b,  where  b  is  the 
impact  parameter.  Figure  6  shows  the  opacity  function  for  the 
three  different  reaction  channels  at  10  km  s_1  using  the  AM1- 
SRP1  surface.  The  internal  energy  deposition  is  correlated  with 
the  distribution  of  impact  parameters:  nearly  all  H-elimination 
reactions  (red)  are  hard  collisions  with  impact  parameters  less 
than  1.0  A.  Nearly  all  CH3  elimination  reactions  occur  for 
impact  parameters  less  than  2  A,  and  the  H-atom  abstraction 
pathway  has  a  broad  range  of  impact  parameters,  mostly 
weighted  between  2  and  3  A. 


Figure  6.  Opacity  function  divided  by  the  bin  width,  P(b)/db,  in  units 
of  (l/A)  calculated  at  a  collision  velocity  of  10  km  s-1  using  the  AM1- 
SRP1  potential.  The  hydrogen  abstraction  reaction  is  shown  in  black, 
the  hydrogen  elimination  reaction  is  shown  in  red,  and  the  methyl 
elimination  reaction  is  shown  in  blue.  Error  bars  indicate  standard 
deviation,  and  theimpact  parameter,  b,  bin  size  is  0.2  A. 


Nonreactive  collisional  excitation  of  DMMP  leads  to  cross 
sections  much  larger  than  reactive  pathways  and  may  carry 
important  collision  dynamics  information.  For  the  nonreactive 
collision  process,  flux— velocity  maps  provide  insight  into  the 
energy  transfer  process.  Figure  7  shows  flux— velocity  contour 
maps  for  nonreactive  collisions  of  0(3P)  +  DMMP  at  velocities 
of  4,  6,  8,  and  10  km  s_1  using  the  AM1-SRP1  Hamiltonian.  For 
reference,  inner  and  outer  circles  are  shown  at  5  and  10  km  s_1, 
respectively.  Forward  scattering  (zero-degree  scattering  angle) 
points  to  the  right,  and  the  contour  color  is  on  a  log  scale  to 
bring  out  important  features.  From  these  features,  there  are  two 
main  scattering  mechanisms,  a  dominant  forward  scattering 
component  with  little  internal  excitation  and  a  lesser  but 
significant  side  and  backscatter  component  which  leads  to 
internal  excitation.  In  addition,  there  is  a  third  mechanism  at 
low  scattering  angles  at  8  and  10  km  s_1,  which  corresponds  to 
large  DMMP  internal  energies.  This  may  correspond  to  O- 
atoms  undergoing  low-impact  parameter  collisions  that  get 
funneled  into  the  forward  scattering  direction. 

Vibrational  Cross  Sections.  For  H-atom  abstraction,  the 
OH  vibrational  state  cross  section  is  plotted  in  Figure  8  at  the 
AMI,  SRP1,  and  SRP2  levels  at  collision  velocities  of  6,  8,  and 
10  km  s_1.  The  abstraction  reactions  from  the  AMI  calculations 
result  in  the  most  vibrationally  excited  OH,  which  is  consistent 
with  the  overly  large  exothermicities  at  this  level  of  theory.  For 
SRP1  and  SRP2,  the  OH  vibrational  distributions  peak  at  n  =  1. 
These  vibrational  distributions  remain  fairly  constant  over  a 
broad  range  of  collision  velocities. 

For  H-abstraction,  the  DMMP  radical  product  internal 
energy  vibrational  and  rotational  cross  sections  at  6,  8,  and  10 
km  s_1  are  shown  in  Figure  9  at  the  AMI,  SRP1,  and  SRP2 
levels  of  theory.  Zero-point  energy  is  subtracted  from  the 
vibrational  binning  in  order  to  isolate  product  internal  energy 
change.  We  have  also  plotted  DMMP  radical  vibrational 
thermal  distributions,  scaled  by  the  total  reaction  cross  section. 
For  6  km  s_1,  the  distributions  plotted  are  800  K  (AMl),  900  K 
(AMl-SRPl),  and  800  IC  (AM1-SRP2);  for  8  km  s"1,  the 
distributions  plotted  are  900  IC  (AMl),  1000  IC  (AMl-SRPl), 
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Figure  7.  Velocity— flux  contours  log10  dc(d£2  dv)  1  in  units  of  cm2  (sr  cm  s  1 )  1  for  the  differential  cross  section  of  nonreactive  collisions  of  DMMP 
+  0(3P)  at  (a)  4,  (b)  6,  (c)  8,  and  (d)  10  km  s-1.  The  results  are  shown  for  the  AM1-SRP1  potential. 


and  800  K  (AM1-SRP2);  and  for  10  km  s  l,  the  distributions 
plotted  are  900  K  (AMI),  900  K  (AMl-SRPl),  and  900  K 
(AM1-SRP2).  These  thermal  distributions  fit  the  calculated 
distributions  fairly  well,  except  for  indications  of  a  high-energy 
tail  at  the  higher  velocities.  For  the  AMl-SRPl,  the  vibrational 
distribution  is  shifted  to  a  higher  energy  which  is  particularly 
noticeable  at  a  collision  velocity  of  8  km  s-1  and  is  consistent 
with  the  average  internal  energy  values  in  Table  3. 

The  vibrational  cross  sections  of  the  DMMP— O  radical 
following  hydrogen  elimination  at  8  and  10  km  s-1  are  shown 
in  Figure  10.  Zero-point  energy  is  subtracted  from  the 
vibrational  binning  in  order  to  isolate  product  internal  energy 
change.  Thermal  vibrational  distributions  (scaled  by  the  cross 
section)  are  plotted  with  dashed  curves.  At  8  km  s-1, 
distributions  of  1700  K  (AMl),  1600  K  (AMl-SRPl),  and 
1700  K  (AM1-SRP2)  are  plotted.  The  thermal  distributions 
plotted  for  10  km  s"1  are  2200  IC  (AMl),  2200  IC  (AMl-SRPl), 
and  2200  K  (AM  1-SRP2).  The  vibrational  cross  sections  for  the 
DMMP— O  radical  product  are  well  fit  by  the  thermal 
distributions,  which  correspond  to  extremely  high  temper¬ 
atures,  and  strongly  depend  on  the  collision  energy.  These  high 
temperature  fits  are  consistent  with  the  large  average  internal 
energy  values  of  Table  4.  Vibrational  cross  sections  for  methyl 
elimination  products  were  not  examined  due  to  the  insufficient 
number  of  methyl  elimination  reactions  to  generate  significant 
distribution  statistics. 

To  characterize  the  nonreactive  distributions,  in  Figure  11  we 
show  the  vibrational  cross  sections  versus  DMMP  internal 
energy  at  4,  8,  and  10  km  s-1  for  the  nonreactive  collision  case, 
where  we  have  subtracted  the  initial  DMMP  zero-point  energy. 


Following  the  analysis  of  the  flux  velocity  maps  of  Figure  7, 
where  forward  and  backward  scattering  mechanisms  were 
found,  we  decompose  the  cross  sections  into  angular  scattering 
components.  At  4  km  s-1,  the  distribution  is  divided  into 
scattering  between  5°  and  70°  and  scattering  greater  than  70°. 
At  8  and  10  km  s-1,  the  distribution  is  divided  into  three  parts: 
5—20°,  20—70°,  and  70—180°.  For  reference,  thermal  vibra¬ 
tional  distributions  of  DMMP  are  shown  in  the  plots  with 
dashed  lines  and  are  scaled  by  the  cross  section  at  a  particular 
range  of  scattering  angles.  For  all  velocities  the  large-angle 
scattering  component  (70—180°)  agrees  well  with  thermal 
distributions  at  400,  1200,  and  1600  K  for  the  4,  8,  and  10  km 
s"1  collision  velocities,  respectively.  At  8  and  10  km  s-1,  the 
agreement  of  the  20—70°  component  with  lower  temperature 
thermal  distributions  is  not  as  good,  with  a  high-energy  tail  in 
the  calculated  results  that  is  not  accounted  for. 

Rotational  Cross  Sections  of  Reactive  Collisions.  For 
the  H-atom  abstraction  channel,  the  rotational  state  cross 
sections  for  OH  at  velocities  of  6,  8,  and  10  km  s-1  are  plotted 
in  Figure  12a,  b,  and  c,  respectively.  The  average  rotational 
energy  increases  and  the  rotational  distributions  widen 
modestly  with  collision  velocity.  Unlike  the  vibrational  cross 
sections,  the  rotational  cross  sections  did  not  fit  Boltzmann 
distributions  well.  Therefore,  to  characterize  the  rotational 
distributions,  we  fit  the  rotational  energy  cross  sections  to 
distributions  derived  from  surprisal  theory 

c(Erot)  ‘X  °°(£rot)  exP(-0/rot) 

where  the  prior  distribution,  <7°,  is  proportional  to  (2/  +  l)(Etot 
— £rot)1/2,  /rot  is  the  energy  fraction  in  rotation,  and  6  is  the 
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Figure  8.  Vibrational  state  cross  sections  for  the  OH  product  using 
AMI  (black),  AM1-SRP1  (red),  and  AM1-SRP2  (blue)  at  collision 
velocities  of  (a)  6,  (b)  8,  and  (c)  10  km  s_1.  Error  bars  indicate 
standard  deviation. 


surprisal  parameter.  These  fits  are  shown  in  dashed  curves  in 
Figure  12  for  the  OH  product  of  the  hydrogen  abstraction 
reactions.  For  the  OH  product,  the  surprisal  parameters  are 
similar  for  all  the  potentials. 

For  the  H-abstraction  channel,  the  rotational  cross  sections 
for  the  residual  DMMP  radical  product  are  given  in  Figure 
13a— c  for  collision  velocities  of  6,  8,  and  10  km  s-1.  The  AMI 
rotational  distributions  are  generally  more  excited  than  the  SRP 
distributions.  This  is  probably  due  to  errors  in  the  AMI 
potential  in  the  product  channel  leading  to  overly  large  AMI 


Figure  9.  Vibration  cross  sections  as  a  function  of  vibrational  energy  in 
kcal  mol-1  for  the  residual  DMMP  radical  product  following  hydrogen 
abstraction  using  AMI  (black),  AM1-SRP1  (red),  and  AM1-SRP2 
(blue)  at  collision  velocities  of  (a)  6,  (b)  8,  and  (c)  10  km  s-1.  The 
dashed  curves  are  thermal  distributions  and  for  comparison  are  scaled 
by  the  abstraction  reaction  cross  section  for  AMI  (black),  AM1-SRP1 
(red),  and  AM1-SRP2  (blue),  respectively.  For  6  km  s-1,  thermal 
distributions  of  800  IC  (AMl),  900  K  (AMl-SRPl),  and  800  K  (AM1- 
SRP2)  are  shown.  For  8  km  s-1,  thermal  distributions  of  900  K 
(AMI),  1000  IC  (AMl-SRPl),  and  800  IC  (AM1-SRP2)  are  shown. 
For  10  km  s-1,  thermal  distributions  of  900  IC  (AMl),  900  IC  (AMl- 
SRPl),  and  900  IC  (AM1-SRP2)  are  shown.  In  all  plots,  error  bars 
which  indicate  standard  deviation  are  given  for  every  third  data  point, 
and  2  kcal  mol-1  energy  bins  are  used. 

exothermicities.  The  SRP-1  cross  section  at  low  rotational 
energies  is  also  significantly  larger  than  the  SRP-2  cross  section. 
This  reflects  a  qualitative  difference  between  the  two  surfaces 
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Figure  10.  Vibration  cross  sections  as  a  function  of  energy  in  kcal 
mol-1  for  the  residual  DMMP  product  following  hydrogen  elimination 
using  AMI  (black),  AM1-SRP1  (red),  and  AM1-SRP2  (blue)  at 
collision  velocities  of  (a)  8  and  (b)  10  km  s-1  .  The  dashed  curves  are 
thermal  distributions  and  for  comparison  are  scaled  by  the  elimination 
reaction  cross  section  for  the  AMI  (black),  AM  1-SRP1  (red),  and 
AM1-SRP2  (blue),  respectively.  For  8  km  s-1,  thermal  distributions  of 
1700  IC  (AMI),  1600  IC  (AMl-SRPl),  and  1700  IC  (AM1-SRP2)  are 
shown.  For  10  km  s-1,  thermal  distributions  of  2200  K  (AMl),  2200  K 
(AMl-SRPl),  and  2200  K  (AM  1-SRP2)  are  shown.  In  all  plots,  error 
bars  which  indicate  standard  deviation  are  given  for  every  third  data 
point,  and  4  kcal  mol-1  energy  bins  are  used. 

seen  also  in  the  vibrational  distributions  and  internal  energies. 
The  AM1-SRP2  surface  is  a  better  overall  fit  to  the  benchmark 
CBS-QB3  results  for  this  pathway,  as  indicated  in  Table  2,  and 
we  believe  the  AM1-SRP2  cross  sections  are  likely  more 
accurate.  The  surprisal  fits  of  the  rotational  distributions  are 
shown  as  dashed  lines  in  Figure  13.  In  this  case,  the  prior 
distributions  were  obtained  by  modeling  the  DMMP  radical  as 
an  asymmetric  top  and  determining  the  rotational  degeneracy 
as  a  function  of  rotational  energy.  The  surprisal  parameter 
values  for  the  unmodified  AMI  and  AMl-SRPl  are  similar.  The 
surprisal  parameter  values  for  the  AM1-SRP2  fits  are  nearly 
three  times  greater  than  both  the  unmodified  AMI  and  AMl- 
SRPl  parameters.  The  rotational  cross  sections  are  given  for 
the  DMMP— O  radical  product  following  H-elimination  at  8 
and  10  km  s-1  in  Figure  14a  and  b,  respectively.  There  is  a 
broad  distribution  of  rotational  energies,  reflecting  the 
unusually  large  deposition  of  initial  translation  into  internal 
energy.  The  surprisal  fits  plotted  in  Figure  14  as  dashed  lines 


Figure  11.  Vibrational  cross  sections  as  a  function  of  energy  in  kcal 
mol-1  for  DMMP  following  a  nonreactive  collision  at  (a)  4,  (b)  8,  and 
(c)  10  km  s-1  using  AMl-SRPl.  The  distribution  at  4  km  s-1  is 
separated  into  two  distributions,  5—70°  (black)  and  70—180°  (blue), 
while  the  distributions  at  8  and  10  km  s-1  are  separated  into  three 
distributions,  5—20°  (black),  20—70°  (red),  and  70—180°  (blue).  The 
dashed  curves  are  thermal  distributions  scaled  by  the  inelastic  cross 
section  for  the  angular  distributions.  For  4  km  s-1,  thermal 
distributions  of  100  K  (5—70°)  and  400  K  (70—180°)  are  shown. 
For  8  km  s-1,  thermal  distributions  of  200  K  (5—20°),  600  K  (20— 
70°),  and  1200  K  (70—180°)  are  shown.  For  10  km  s-1,  thermal 
distributions  of  300  IC  (5-20°),  800  IC  (20-70°),  and  1600  IC  (70— 
180°)  are  shown.  In  all  plots  error  bars  which  indicate  standard 
deviation  are  given  for  every  third  data  point,  and  2  kcal  mol-1  energy 
bins  are  used. 
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Figure  12.  Rotational  state  cross  sections  as  a  function  of  energy  in 
kcal  mol-1  for  the  OH  product  using  AMI  (black),  AM1-SRP1  (red), 
and  AM1-SRP2  (blue)  at  collision  velocities  of  (a)  6,  (b)  8,  and  (c)  10 
km  s_1.  The  dashed  curves  are  surprisal  fits  for  the  rotational 
distributions:  (a)  =  1.7,  $ami-srpi  =  2.0,  $ami-spr2  =  1.6;  (b) 

-  3.5,  #ami-srpi  -  3.4,  0ami-spr2  =  2.8;  (c)  $ami  -  4.6,  $ami-srpi  =  4.0, 
0ami-spr2  -  4.7.  In  all  plots,  error  bars  which  indicate  standard 
deviation  are  given  for  every  third  data  point,  and  4  kcal  mol-1  energy 
bins  are  used. 


for  each  potential  are  similar  for  all  three  potentials.  As  with  the 
vibrational  distributions,  the  methyl  elimination  rotational 
distribution  statistics  were  too  poor  to  effectively  analyze. 

DMMP  Spectra.  Optical  detection  of  0(3P)  +  DMMP 
nonreactive  collisions  may  be  an  important  probe  in  future 
hyperthermal  beam  measurements.  As  suggested  by  Figures  7 
and  1 1  and  in  analysis  of  the  nonreacted  rotational  distributions 


Figure  13.  Rotational  cross  sections  as  a  function  of  energy  in  kcal 
mol-1  for  the  residual  DMMP  radical  product  following  hydrogen 
abstraction  using  AMI  (black),  AM1-SRP1  (red),  and  AM1-SRP2 
(blue)  at  collision  velocities  of  (a)  6,  (b)  8,  and  (c)  10  km  s-1.  The 
dashed  curves  are  surprisal  fits  for  the  rotational  distributions:  (a)  d^i 
-  7.7,  $ami-srpi  “  8.6,  $ami-spr2  =  25.3;  (b)  $ami  -  13.9,  $ami-srpi  - 
15.9,  $ami-spr2  =  44.4;  (c)  0AMl  =  21.9,  $ami-srpi  =  25.1,  $ami-spr2  = 
72.8.  In  all  plots,  error  bars  which  indicate  standard  deviation  are  given 
for  every  third  data  point,  and  4  kcal  mol-1  energy  bins  are  used. 


which  we  have  not  shown,  product  internal  excitation  will  be 
sensitive  to  the  scattering  angle  and  can  be  characterized  to 
some  extent  by  a  Boltzmann  temperature.  High-temperature 
spectra  of  DMMP  are  not  available.  Therefore,  we  have 
simulated  DMMP  at  several  representative  temperatures  (298, 
800,  and  1600  K)  which  span  the  general  expected  range  of  the 
product  distributions.  These  spectra  are  plotted  in  Figure  15. 
They  were  obtained  by  a  Monte  Carlo  technique  which  uses 
known  dipole  transition  moments  and  other  spectroscopic 
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Figure  14.  Rotational  cross  sections  as  a  function  of  energy  in  kcal 
mol-1  for  the  residual  DMMP  product  following  hydrogen  elimination 
using  AMI  (black),  AM1-SRP1  (red),  and  AM1-SRP2  (blue)  at 
collision  velocities  of  (a)  8  and  (b)  10  km  s_1.  The  dashed  curves  are 
surprisal  fits  for  the  rotational  distributions:  (a)  d^i  ~  4.7,  $ami-srpi  = 
6.8,  $ami-spr2  =  3.2;  (b)  =  2.9,  $ami-srpi  =  3.5,  0ami-spr2  =  3.6.  In 

all  plots,  error  bars  which  indicate  standard  deviation  are  given  for 
every  third  data  point,  and  4  kcal  mol-1  energy  bins  are  used. 


Figure  15.  Calculated  spectra  of  DMMP  with  T  =  298  K  (black),  T  = 
800  K  (red),  and  T  =  1600  K  (blue).  The  spectrum  from  the  NIST 
database  at  373  K  is  also  shown  (green). 

information  to  generate  representative  spectral  lines.37  Large 
numbers  of  individual  ro-vibrational  line  shapes  are  generated, 


typically  tens  of  thousands,  depending  on  temperature.  These 
representative  lines  are  weighted  in  a  Monte  Carlo  fashion,  and 
each  vibrational  band  is  scaled  by  its  dipole  strength.  A 
spectrum  is  built  up  which  yields  a  representation  of  the  true 
spectrum  which  can  contain  millions  of  lines  or  more.  The 
DMMP  373  K  spectrum  from  the  NIST  database38  is  also 
plotted  in  Figure  15  for  comparison.  The  calculated  spectra  are 
normalized  by  the  integrated  intensity  of  the  NIST  spectrum. 
Additionally,  the  vibrational  modes  with  strong  features  are 
indicated  in  the  figure.  As  expected,  the  spectrum  broadens 
with  temperature.  As  indicated  in  the  figure,  a  measurement  of 
the  spectral  peak  widths  in  a  crossed-beam  apparatus  should 
yield  an  estimate  of  the  nonreactive  collision  product 
temperatures. 

IV.  CONCLUSIONS 

We  developed  two  global  potential  energy  surfaces  for  0(3P)  + 
DMMP  collisions,  using  the  specific  reaction  parameters  (SRP) 
approach.  Each  global  surface  is  simultaneously  fit  along  the 
three  major  reaction  pathways:  hydrogen  abstraction,  hydrogen 
elimination,  and  methyl  elimination.  We  used  these  global 
potential  surfaces  in  classical  dynamics  simulations  and 
computed  reactive  cross  sections  from  4  to  10  km  s-1  collision 
velocity.  We  examined  the  energy  disposal  and  angular 
distributions  of  the  reactive  and  nonreactive  collided  products. 
We  found  that  for  reactive  collisions,  an  unusually  large  amount 
of  the  initial  collision  energy  is  transformed  into  internal 
energy.  We  also  found  that  internal  energy  deposition  is 
correlated  with  the  distribution  of  impact  parameters:  nearly  all 
H-elimination  reactions  are  hard  collisions  with  impact 
parameters  less  than  1.0  A,  and  about  80%  of  available  energy 
is  deposited  in  the  internal  energy  of  the  products.  Nearly  all 
CH3  elimination  reactions  occur  for  impact  parameters  less 
than  2  A  and  result  in  about  60—70%  of  available  energy 
turning  into  internal  energy,  and  the  H-atom  abstraction 
pathway  has  a  broad  range  of  impact  parameters  up  to  about  3 
A  and  results  in  about  50—60%  of  available  energy  transformed 
into  internal  product  energy.  For  the  elimination  reactions, 
most  of  this  internal  excitation  goes  into  vibration.  We  also 
analyzed  the  nonreactive  and  reactive  product  internal  energy 
distributions,  many  of  which  fit  Boltzmann  temperatures  up  to 
~2000  K,  and  predict  excited  product  spectra. 
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